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ABSTRACT 

Microlensing in Q2237+0305 between 1985 and 1995 (e.g. Irwin et al. 1989; Cor- 
rigan et al 1991; 0stensen et al. 1996) has been interpreted in two different ways. 
Firstly, the observed variations can be explained through microlensing by stellar mass 
objects of a continuum source with dimensions significantly smaller than the microlens 
Einstein Radius {rjo), but consistent with that expected for thermal accretion discs 
(e.g. Wambsganss, Paczynski & Schneider 1990; Ranch & Blandford 1991). However, 
other studies have shown that models having sources as large as 5 rjo can reproduce the 
observed variation (Refsdal & Stabell 1993; Haugan 1996). In this paper we present 
evidence in favour of a small source. Our approach uses the distribution of microlensed 
light-curve derivatives to place statistical limits (as a function of source size) on the 
number of microlens Einstein radii crossed by the source during the monitoring period. 
In contrast to previous analyses, our results are therefore not dependent on an assumed 
time-scale. Limits on the source size are obtained from two separate light-curve fea- 
tures. Firstly, recently published monitoring data (Wozniak et al. 2000a, b; OGLE web 
page) shows large variations (~.8-1.5 magnitudes) between image brightnesses over a 
period of ~ 700 days or ~15% of the monitoring period. Secondly, the 1988 peak in the 
image A light-curve had a duration that is a small fraction (<0.02) of the monitoring 
period. Such rapid microlensing rises and short microlensing peaks only occur for small 
sources. We find that the observed large-rapid variation limits the source size to be 
<Q.2rjo (95% confidence). The width of the light-curve peak provides a stronger con- 
straint of <0. 025770 (99% confidence). The Einstein radius (projected into the source 
plane) of the average microlens mass (m) in Q2237+0305 is rjo ^ 10^^-\/(m) cto. The 
interpretation that stars are responsible for microlensing in Q2237+0305 therefore re- 
sults in limits on the continuum source size that are consistent with current accretion 
disc theory. 

Key words: gravitational lensing - microlensing - numerical methods. 



1 INTRODUCTION 

Tlie object Q2237+0305 (Huchra et al. 1985) comprises a 
source quasar at a redshift of 2: = 1.695 that is gravitation- 
ally lensed by a foreground galaxy with z = 0.0394 produc- 
ing 4 resolvable images with separations of ~ 1". Each of the 
4 images are observed through the galactic bulge, which has 
a microlensing optical depth in stars that is of order unity 
(e.g. Kent & Falco 1988; Schneider et al. 1988; Schmidt, 
Webster & Lewis 1998). In addition, the proximity of the 
lensing galaxy means that the effective transverse velocity 
may be high, yielding an expected microlensing event time- 
scale significantly shorter than that of other lensed quasars. 
The combination of these considerations make Q2237-f 0305 
the ideal object from which to study microlensing. Indeed, 



Q2237+0305 is the only object in which cosmological mi- 
crolensing has been directly confirmed (Irwin et.al 1989; 
Corrigan et.al 1991; Wozniak et al. 2000a,b). 

Initially, this confirmation came in the form of a ~0.2 
magnitude brightening of image A with a rise-time of ~ 26 
days (Corrigan et al. 1991). Wambsganss, Paczynski & 
Schneider (1990) found that, assuming a galactic transverse 
velocity of ~ 600fcmsec~^, this rise could be explained by 
microlensing due to stellar masses of a source having dimen- 
sions much (< 0.01 ER) smaller than the microlens Einstein 
radius (ER) and therefore the typical caustic spacing. 

Spectral observations of subsequent microlensing also 
support the case for a source that is small with respect to 
the microlens ER. Lewis et al. (1998) determined the ra- 
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tios of emission line equivalent widths relative to one image. 
They show that the ratios vary between images at a sin- 
gle epoch, and that the ratio for a single image (image A) 
varies between two different epochs. This spectral change is 
interpreted as being due to the different spatial extents of 
the continuum and emission line regions being differentially 
amplified due to microlensing, and suggests that the contin- 
uum region is smaller than the typical caustic separation. 

Refsdal & Stabell (f993) proposed a model with a very 
low average microlens mass ('^ 1Q~'^ Mq), and a source size 
consistent in physical size with the models of Wambsganss, 
Paczynski & Schneider (f990). In this model the source size 
is several ER. Surprisingly, brightness variations as large 
as ~ 0.5 magnitudes are predicted by this model. On the 
other hand Witt & Mao (1994) note that smeared out light 
curves are produced which have trouble producing asym- 
metric events such as that observed for image A (1989-90), 
although this interpretation was disputed by Haugan (1996). 
Very large sources, (> 7ER) have been ruled out by Refsdal 
and Stabell (1997) from the CIV line data of Lewis et al 
(1998). 

There are two important reasons to distinguish between 
the cases of a source that is large/small with respect to the 
microlens ER. Firstly, much observational effort (OGLE col- 
laboration; Lens Monitoring Project, Apache Point Obser- 
vatory) is currently being directed towards monitoring of 
Q2237+0305 in the hope of observing a caustic crossing (the 
event time-scale is of order months rather than years due to 
the large distance ratio). As has been discussed by several 
authors (e.g. Grieger et al. 1988; Grieger et al. 1991; Agol & 
Krolic 1999; Mineshige & Yonehara 1999), the light curve of 
a straight-single caustic event contains information on the 
source geometry on nano-arcsecond scales, providing oth- 
erwise unobtainable resolution for observation of a quasars 
continuum region. The mode of analysis discussed by these 
authors is only valid in the case of a source that is much 
smaller than the typical caustic spacing. In addition, the 
differential magnification probed by these methods is much 
more significant for small sources. Secondly, a good under- 
standing of the microlensing parameters (e.g. galactic trans- 
verse velocity and mean microlens mass) are required for the 
successful analysis of a caustic crossing event. One means of 
obtaining such understanding is to interpret the monitoring 
data in terms of the microlensing rate (e.g. Lewis & Irwin 
1996; Wyithe, Webster & Turner 1999, 2000b, (hereafter 
WWT99, WWTOOb)). The rate is approximately indepen- 
dent of source size for small sources, a feature that can be 
used to remove one degree of freedom from the problem. 

This paper presents arguments that support the hy- 
pothesis of a source that has dimensions smaller than the 
ER. Sections ^ and ^ describe the microlensing models and 
the collection of published monitoring data, while section 
describes how the distribution of microlensed light-curve 
derivatives can be used to place limits on the length of caus- 
tic network sampled by the observations. In section ^ we 
discuss how two different light-curve features, in combina- 
tion with the sampling length, limit the continuum source to 
be significantly smaller than the microlens ER. This paper 
differs from previous analyses by concentrating on determin- 
ing the ratio of the source size to the microlens mass, rather 
than the source size in physical units. A more important dif- 
ference is that, unlike previous work, our analysis contains 



Table 1. Values of the total optical depth (k) and the magni- 
tude of the shear (7) at the position of each of the 4 images of 
Q2237+0305. The quoted values are from Schmidt, Webster & 
Lewis (1998). 



Image 




I7I 


A 


0.36 


0.40 


B 


0.36 


0.40 


C 


0.69 


0.71 


D 


0.59 


0.61 



no assumptions about the mean microlens mass or the trans- 
verse velocity (other than a prior probability for transverse 
velocity which, within reasonable limits, has no effect on our 
resuhs (WWTOOb)). 



2 THE MICROLENSING MODEL 

To model microlensing in Q2237-f0305 we assume the 
macro-parameters for the lensing galaxy calculated by 
Schmidt, Webster & Lewis (1998). These values are shown 
in table |l|. Two orientations were chosen for the transverse 
velocity with respect to the galaxy, with the source trajec- 
tory being parallel to the A— B and C— D axes. Photometric 
errors were simulated by applying a random perturbation to 
the model light-curves having a magnitude distributed ac- 
cording to a Gaussian of halfwidth a. The simulations used 
two different estimates of the error in the photometric mag- 
nitudes. In the first case a small error was assumed (SE). 
For images A and B, (ts_b=0.01 mag, and for images C and 
D asE=0.02 mag. In the second case, a larger error was as- 
sumed (LE). For images A and B, cr_t_B=0.02 mag and for 
images C and D (j_l_e=0.04 mag. The observational error 
in Irwin et al. (1989) was 0.02 mag. These models produce 
quantitatively similar results. Therefore, in this paper we 
present only results from models with small errors (SE) and 
a transverse velocity whose direction lies along the A-B axis. 

Both the microlensing rate due to a transverse velocity 
(e.g. Witt, Kaiser & Refsdal 1993; Lewis & Irwin 1996), as 
well as the corresponding rate due to random proper motions 
(Wyithe, Webster & Turner 2000a (hereafter WWTOOa)) are 
not functions of the details of the microlens mass distribu- 
tion, but rather are dependent only on the mean microlens 
mass. The independence of the microlensing rate on the form 
of the mass-function has been checked for models contain- 
ing mass ranges over 2 orders of magnitude. We assume that 
the dominant contribution to the optical depth comes from 
objects with masses differing by less than 2 orders of mag- 
nitude and consequently limit our attention to models in 
which all the point masses have the same mass since the 
results obtained will be applicable to other models with dif- 
ferent forms for the mass function. In addition, we consider 
only models that contain no continuously distributed mat- 
ter. 

The distribution of microlensed light-curve derivatives 
has nearly the same form (up to a scaling factor in the 
derivative) whether the variation results from a transverse 
velocity of a static screen of stars, or whether it is due to 
their random proper motion (WWTOOa). This suggests that 
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microlensing due to the combination of a transverse veloc- 
ity and stellar proper motions can be approximated using a 
static screen of microlenses with an effective transverse ve- 
locity (which is larger than the physical galactic transverse 
velocity), thus the microlensing effect of microlens proper 
motions is included implicitly (WWT99). In this paper we 
use model light-curves obtained using this effective trans- 
verse velocity approximation. 

The ER of a microlens in the source plane is denoted 
by rjo- We consider Gaussian sources with sizes S (where S 
is the half width of the Gaussian intensity profile) ranging 
over three orders of magnitude, from S — O.OOISt^o — 1.6?;o- 
To construct microlensing light-curves for S < O.OSiyo we 
use the inversion technique of Lewis et al. (1993) and Witt 
(1993). For S > O.lrjo we use the ray-tracing method (e.g. 
Kayser, Refsdal & Stabell 1986; Schneider & Weiss 1987; 
Wambsganss, Paczynski & Katz 1989). 

2.1 Models for small sources 

The microlensing models for sources of size S < O.OSryo 
presented in this work have been discussed in detail in 
WWTOOb. Finite source light-curves were produced by con- 
volving a Gaussian source profile with a point-source light- 
curve (e.g. Witt & Mao 1994). For each combination of mi- 
crolensing parameters, 100 light-curves of length lOrjo were 
produced. 

2.2 Models for large sources 

For source sizes S > O.lryo, the 1-D approximation to the 
finite source light-curve is no longer valid, and the lower 
resolution requirements make ray-tracing the appropriate 
method. For sources having OArjo < S < OArjo, light-curves 
were produced from magnification maps having side-lengths 
of AOrjo. For sources with sizes 0.8r;o < S < l-Grjo, we 
produced magnification maps having side- lengths of 160r;o. 
The magnification maps had a resolution of 500x500 pix- 
els, and the number of stars used in the models was calcu- 
lated through the method described in Lewis & Irwin (1995) 
and Wyithe & Webster (1999). Finite source light-curves are 
produced from these maps through convolution with a 2-d 
Gaussian source profile. 



3 MONITORING DATA 

We have compiled a data set that includes the photometry 
presented in Schneider et al. (1988), Kent & Falco (1989), 
Irwin et al. (1989), Corrigan et al. (1991), 0stensen et al. 
(1996); Wozniak et al. (2000a,b) and the OGLE web page 

(see 

http: / / www.astro.princeton.edu/~ogle/ogle2 /huchra.htm] ) . 
To successfully analyse the data, even sampling and accurate 
photometry are required so that light-curve derivatives at 
one or a few epochs do not dominate the statistics. We have 
therefore averaged points on the light-curves, producing a 
more even sampling rate and reducing photometric error ac- 
cording to the procedure described in WWT99. Points were 
averaged if the observations were taken within one week, 
a procedure that does not smooth out observed short fluc- 
tuations (less than 1 week) since these are not present in 



the data above the noise generated by photometric uncer- 
tainty. Following this, any points having an associated error 
above AM = 0.05 magnitudes were removed from the sam- 
ple since data points with large errors substantially degrade 
the measurement of microlensing rate through introduction 
of noise into the low derivative regime. There were also two 
data points that were discarded since they displayed corre- 
lated flux variation in two images. The upper plot in figure 
^ shows the entire data set as well as the collection of points 
used for our analysis, dotted sections highlight large gaps be- 
tween observations. Error bars are shown representing the 
published errors. Errors for averaged points have been added 
in quadrature. 

While data taken prior to ~ 3100 days is in R and r- 
bands, observations from OGLE have been made in V-band. 
We have not included a colour correction since our analysis 
considers the difference between image magnitudes. In addi- 
tion, the continuum at different wavelengths may be emit- 
ted from regions having different physical scales. Differen- 
tial magnification due to microlensing therefore introduces a 
complication when comparing data taken in different bands. 
However, microlens induced colour change occurs primarily 
during caustic crossings (Wambsganss & Paczynski 1991). 
Therefore if S ^ 770, colour change is only observed for a 
small fraction of time, while if S ~ 770 the differences are 
small. The lower plot of figure ^ shows the set of difference 
light-curves on which our analysis is performed. We have 
drawn solid lines between all points on this plot since the 
derivatives between all points are used in the analysis. 



4 THE SAMPLING LENGTH 

WWT99,WWT00b describe a method to compute the prob- 
ability for the value of the lensing galaxies effective trans- 
verse velocity using the distribution of microlensed light- 
curve derivatives. The effective transverse velocity obtained 
is in uints of kmsec"^ {m)~2 where (m) is the mean mi- 
crolens mass. When multiplied by the length of the monitor- 
ing period, this yields the length of caustic structure sam- 
pled ijiperiod) in units of the ER of the average microlens 
mass (?7o). This is determined as a function of source size 
S. Note that r]period{S) is free of assumptions about time- 
scale and microlens mass. We have computed the probabil- 
ity Pr^{rjperiod\S) for sourcc slzes 5* = 0.0015625 0.003125, 
0.00625 0.0125, 0.025, 0.05, 0.1, 0.2, 0.4, 0.8 and l.Qr)o. For 

5 < O.lrjo and 5* > O.lrjo, Pr]{'riperiod\S) was computed from 
5000 and 1000 mock sets of monitoring data respectively 
at each of 50 assumed sampling lengths rj. The cumulative 
distributions 

PviV < Vperiod) ^ / Pv('n'\S)dT]' (1) 

Jo 

are shown in the upper plot of figure ^ (thicker lines in- 
dicate larger sources). For S < O.O5770 (solid lines), the 
measurement of sample length is approximately independent 
of source size (this independence is discussed in WWT99). 
However, the measured sample length increases as larger 
source sizes are assumed because larger sources have de- 
creased event amplitudes and increased event time-scales. 
This combination results in a higher transverse speed being 
required to reproduce the observed derivatives. The limits 
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Figure 1. Compilation of published light-curve data for Q2237+0305. Top: Individual image light-curves. The light lines show the entire 
data set, and the dark lines show the modified data-set used for our analysis. Dotted lines span large gaps in the observations. Bottom: 
Six difference light-curves (A-B, A-C, A-D, B-C, B-D, C-D) calculated from the modified data set. 



on sample length are ~ 5r;o for small sources, and ~ lOOr^o 
for all sizes considered. 



5 LIMITS ON SOURCE-SIZE 

In this section we place limits on the source size in units 
of microlens ER by combining sampling length limits with 
features in the monitoring light-curves. There are at least 
two features in the monitoring data that suggest a small 
source size. Firstly, while the 0.2 magnitude variation in 
image A (1989) can be explained by a large source, the 



associated peak has a very short duration with respect 
to the monitoring period (~ 2%). The sample length 
iVperiodiS)) must therefore be greater than 50 5* regard- 
less of whether the peak is due to the source passing 
outside of a cusp or to a caustic crossing (the minimum 
event width is set by the source size in either case). We 
show below that a sample length of 50 S can only be 
reconciled with a small source. Secondly, new monitor- 
ing data (Wozniak et al. 2000a,b; OGLE web page (se e 



http:/ /www. astro. princeton.edu/~ogle/ogle2/huchra. html)) 
shows large scale brightness variation in all images, with val- 
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Sample Length (Einstein Radii) 

Figure 2. Top: The probability Pri{ri < r^period I 'S') that the source has crossed less than rip^^-iad ER during the monitoring period. 
Centre: The cumulative probability Pw{W < Wo\S) for light-curve peak width. Bottom: The cumulative probability Pm({AAf}) > 
{ AMoi,s}|S', »7perio£j) for variation in a mock data set being larger than the observed variation. The probabilities arc shown for S = 
0.0015625 0.003125, 0.00625 0.0125, 0.025, 0.05, 0.1, 0.2, 0.4, 0.8, 1.6 rjo (Thicker lines denote larger sources). The solid and dotted lines 
correspond to sources with S < 0.05»yo and S > 0.1»7o respectively. 
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ues ranging from ^.5-1.25 magnitudes over a period of ~700 
days. While we find that (surprisingly) such large microlens- 
ing variation can be produced by large sources over long 
time-scales, rapid large scale microlensing variation only oc- 
curs for small sources. The quantitative analysis of the limits 
imposed on the source size by these features are discussed 
in turn. 

We note that choosing light-curve features a posteriori, 
and comparing them statistically to a sample of models in 
order to draw conclusions regarding input parameters (such 
as the source size) has the potential to introduce a bias. How- 
ever in the present case we feel that our method is justified 
because of the well established prior knowledge that smaller 
microlensed sources produce larger, more rapid variations. 

5.1 Source-size limits from the short peak 

The 1988 peak in the light-curve of image A has been in- 
terpreted as one half of a double horned event (e.g. Racine 
1992), though the data are not conclusive on this point. It 
had a height J> 0.2 magnitudes (the points measure the lower 
limit), and a duration of ^ 100 days (the available points 
measure an upper limit). The duration of the peak is ^ 2.1% 
of the total monitoring period (~ 4700 days). We calculate 
the likely-hood of this duration for different source sizes. 

Firstly, we search for peaks in model light-curves (we 
do not demand that a peak be due to a caustic crossing) 
with maxima more than 0.2 magnitudes above the associ- 
ated minima on both sides. Our peak statistics were ob- 
tained from sampling at the resolution of the simulation 
(>1 point per day). The sample of peaks for sources with 
5* < O.OSt^o was computed from IOOOt^o of light-curve. For 
5* > O.IOt^o, a sample of 1000 light-curve peaks was obtained. 
The central plot of figure |^ shows the cumulative distribu- 
tions Pw(W < Wo\S) of peak fuU-width-at-half-maximums 
(W) given source-size S (sizes assumed are the same as those 
in section^. W is always larger than S and may be 10-100 
times larger. For smaller sources the distribution of W has 
two components. The typical W for a peak resulting from 
a caustic crossing is approximately proportional to S. How- 
ever, if the peak is the result of the source having passed 
outside a cusp then the ratio W : S may be arbitrarily 
large, since in this case W is approximately independent of 
S for 5* ^ 0.05j7o. For larger sources, W's for caustic crossing 
events are typically as large as those for cusp related events, 
while for S ^ r/o no classification is possible since the source 
is generally in contact with two or more caustics. 

From model light-curves we calculate the ratios of W 
to Tjpgriod- The probability P of finding a ratio 



R{S) 



W{S) 



(2) 



^period (5") 

smaller than that observed {Robs) as a function of source 
size is 



P{R<Robs\S) 



p,n(W'\S)pr,{%^riod\S)dr]'p^„oddW' ,{?,) 



''pe,-iad R„bs 



where pw — dPw/dW. Wmax and the largest val- 

ues for W and r]period considered in the simulations and have 
associated probabilities of 0. The left-hand plot of figure ^ 



shows the probability P{R < RobB\S). We find that the short 
duration of the light-curve peak relative to the monitoring 
period limits the source size S to be smaller than ~ 0.025r;o 
at the 99% level. 



5.2 Source size limits from large, rapid 
microlensing variability 

In this section we consider observed large scale variation over 
a period of ~2 years (1997-1999). On this time-scale, intrin- 
sic source variation may be an important factor in the abso- 
lute brightness variation of a single image. We therefore con- 
sider the relative change between image magnitudes, which 
due to the short time-delay (of order 1 day (e.g. Schnei- 
der et al. 1988)) is approximately independent of intrinsic 
variation. The maximum magnitude changes in difference 
light-curves over periods < 700 days (~ 15% of the sample 
length) are AMA-i3=0.83, AMa-c=0.92, AMa-d=1.00, 
AMs_c=1.50, AAfs_i5=0.91 and AMc-i3=0.96. We de- 
fine a set of observed magnitude variations in ascending or- 
der: {AMobs} = {0.83,0.91,0.92,0.96,1.00,1.50}. For each 
assumed monitoring period, 5000 (for 5* < 0.05j7o) or 1000 
(for 5* > O.IO770) mock data sets were calculated using a sam- 
pling rate with identical relative spacing to the monitoring 
data. For each simulated data set we calculate {AAf } such 
that AAf' > AA/'"^ for 2 < i < 6, thus we do not specify 
between which pairs of images the observed changes should 
be seen. We calculate the probability 



Pm 

({AA/} > {/S.Mobs}\S,-qp„ 



(4) 



of finding a change in the difference light-curves of AAf' > 
AAf'ftg magnitudes for all i over O.lSryperiod (during any part 
of the period). The distributions are shown in the lower 
panel of figure ^ Larger sources are less likely to exhibit 
the observed variation over ^15% of a given sample length. 
For example, while model source sizes of 5* O.O5770 always 
exhibit the observed level of microlensing variation over sam- 
ple lengths ^ lOryo, the 5* = l-Grjo source attains these values 
less than 3% of the time over a sample length of ~ 150?7o. 

Pm({AA//} > {A.Moba}\S,-qperiod) IS COUVOlvcd wlth 
Pr]ij)period\S) tO find 

P({AAf} > {AAf^iJIS-) = 



Pm({AA4"} > {/^Mobs}\S,r)period) Prii-qperiodlS) driperijfi) 



This function is shown in the right-hand panel of figure^ We 
find that the large, rapid changes in the image magnitudes 
limit the source to be smaller than O.2770 at the 95% level. 

We note that the calculation of the probabilities 
P({AAf} > {AMobs}\S) assumes that P„({AAf} > 
{AMobs}\S,'npe„od) a.nd Pn{i]period\S) are independent. This 
is a false assumption since the two functions have been 
calculated from the same data set. However, the observa- 
tion of the large scale variation (in the OGLE data) in- 
troduces large derivatives which increases the estimate of 
sample length (over the estimate that would be made from 
the pre 1996 data alone, scaled by the relative monitoring 
periods). In the case of a longer sample length, a source 
has more caustic network over which to undergo large-scale 
microlensing variation. This results in the inference of a 
larger upper limit on source size than would be expected 
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0.01 I 0.1 1 ^ 0.01 0.1 I 1 

Radius Radius 

Figure 3. Left: The probability P{R < Robsl^) of observing a peak width to sample-length ratio smaller than that observed (Robs) a-s 
a function of source size. The dashed line shows the 99% confidence limit. Right: The probability P{{AM} > { AMobsHS) of observing 
changes in the difference light-curves of AM* > AA-P^^ magnitudes for all i during a period of O.lSrjperio^; as a function of source-size. 
The dashed line shows the 95% confidence limit. 



if the sampling length were deduced from the pre-1996 
data alone. Any bias introduced by the co-dependence of 
Pm{{AM} > {AMobs}\S,r]period) and p,,(7?pe„od|5') there- 
fore results in more conservative (ie. larger) source-size lim- 
its than the case where the sample length and large-scale 
variation were measured from different data-sets. 



6 CONCLUSION 

We have used the distribution of microlensed light-curve 
derivatives to find probability functions for the length of 
caustic structure sampled by monitoring observations of 
Q2237-I-0305 as a function of assumed source size. The 1988 
light-curve peak had a height larger than 0.2 magnitudes, 
and a duration less than 0.02 of the monitoring period. At 
the 99% level such a short peak can only be explained for 
a source with dimensions smaller than ~0.025r7o- In addi- 
tion, monitoring by OGLE shows changes in difference light- 
curves ranging from ~ 0.8 — 1.5 magnitudes over ^^15% of 
the monitoring period. We find that such rapid large scale 
changes can only be explained by a source that is smaller 
than ~ 0.2j7o (95% confidence). Importantly, these limits are 
independent of any assumption about mean microlens mass 
or galactic transverse velocity. 

The Einstein radius of the average microlens mass (m) 
in Q2237-f0305 is rjo ~ 10"'"''^ (m) cm. In combination with 
our limits in S, the assumption of stellar mass microlenses 
therefore impose a limit of ^ 2 X lO^'' - 2 x 10^^ cm on the 
continuum source size (consistent with the typical scale-size 
expected for a continuum emitting accretion disc about a 
super-massive black hole (e.g. Rees 1984; Ranch & Bland- 
ford 1991; Agol & Krolik 1999). Conversely, if a source size 
of ~ lO^^cm is assumed, then the our limits on microlens 
mass corresponding to the short light-curve peak and large 
microlens variation are respectively (m) > 0.25AfQ (99%) 
and (m) > O.OO25M0 (95%). Note that for a source having 
a diameter of ^ 10^^ cm we have explicitly checked masses 
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down to (m) ~ 10 ^Mq. A much smaller source must there- 
fore be assumed if sub-stellar masses were to form the bulk 
of the galactic bulge in Q2237-I-0305. 

Our calculations assume that the dominant contribu- 
tion to optical depth comes from a population of microlenses 
having a range smaller than two orders of magnitude. We 
note that our model therefore does not account for the pos- 
sibility suggested by Refsdal & Stabell (1993) in which a 
population of stars and another of very low-mass objects 
both contribute significantly to the optical depth. In this 
case, the rapid variation due to the small masses, and the 
slower variation due to the large masses should be approx- 
imately superimposed (Refsdal & Stabell 1993). However 
our results can be qualitatively interpreted for this scenario. 
There are two possibilities. Firstly, if the source has dimen- 
sions 5* <^ rjsrn (where rjsm is the Einstein radius of the mean 
microlens mass of the low mass population) , then our results 
hold since (m) is quite insensitive to the value of the heavy 
masses. Secondly, if 5 3> rjsm then the light-curves will be 
like those resulting from microlensing by the large masses 
only (with a continuous component of mass density), com- 
bined with a continuous low amplitude flicker. In this case, 
the results we have presented can be interpreted as the ratio 
of source size to the Einstein radius of the mean microlens 
mass of the large population. We therefore conclude that 
a significant contribution to the optical depth comes from 
objects with Einstein radius larger than the source size, re- 
gardless of the mass function. 

The conclusion that the source is small with respect 
to the microlens Einstein radius is important because it 
means that caustic crossing light-curves can be inverted to 
obtain extremely high resolution information on the contin- 
uum source structure. 
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